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We use radial basis functions to model the input-output 
response of an electronic device. A new methodology for pro- 
ducing models that accuratly describe the response of the 
device over a wide range of operating points is introduced. 
A key to the success of the method is the ability to find a 
polynomial relationship between the model parameters and 
the operating points of the device. 

05.45.+b 



I. INTRODUCTION 

In this letter we investigate the possibility of apply- 
ing methods recently developed for modeling nonlinear 
chaotic systems to modelling input-output dynamics of 
a nonlinear electronic device. The ability to describe the 
behaviour of modern electronic devices is becoming in- 
creasingly important. As technology advances, and the 
demands on current technology increase, it is becoming 
more and more apparent that traditional linear modelling 
methods fail to describe the behaviour of modern elec- 
tronic devices . This is particularly true when the de- 
vice is being operated outside of its "normal" (linear) 
response regime, and is subject to complicated (nonlin- 
ear) driving signals [|]||. 

It has been emphasized by researchers in nonlinear dy- 
namics that complicated behaviour, which appears ran- 
dom, can sometimes arise from a low dimensional deter- 
ministic nonlinear system |^|-[|. It seems apt, then, to 
investigate whether nonlinear modelling techniques can 
supplement traditional linear modelling techniques when 
they fail to describe a devices behaviour j7j-||] . 

One result of our study is a new methodology for mod- 
elling electronic devices. In particular, we focus on con- 
structing transportable behavioural models of such de- 
vices. In this research the dynamics of the device is rep- 
resented by input-output data obtained at a few operat- 
ing regimes. A behavioural model is constructed by ad- 
justing its parameters until it accurately mimics the dy- 
namics of the device. Within this context, a behavioural 



model is transportable if it accurately mimics the dynam- 
ics of the device in operating regimes which are different 
from the ones used to construct the model. The ultimate 
goal is a single behavioural model that accurately mimics 
the device in all of its operating regimes. 

Ideally we should compare the performance of our 
models to that of "standard" models using data obtained 
from a real device. However, testing with a real device 
can be expensive. Therefore, in the absence of real data, 
we shall regard the "standard" model for such a device 
as reality. If our behavioural models match the "stan- 
dard" model to within a few percent then we claim that 
we have shown the potential of our methods. 

To illustrate the method we will study a device with 
mild nonlinearities: a simplification of the "standard" 
Ebers-Moll model for a bi-polar junction transistor, 
which we embed in a self-biasing circuit (see Figs. ^ 
and [l]). That is, we are using the transistor as a simple 
amplifier. We will construct transportable behavioural 
models using time domain input-output data obtained 
from the Ebers-Moll model 

In Fig. |^, Vbb and V cc are potential sources. If we ap- 
ply dc-potentials at these points (so called biasing) then 
we are selecting an operating point for the device. A 
graphical summary of this part of the devices behaviour 
is given by the dc V — I characteristics obtained by bias- 
ing the device at different potentials. Figure || shows the 
V — I characteristics of our Ebers-Moll transistor model. 
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FIG. 1. Modified Ebers-Moll transistor model. 
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FIG. 2. Self-biasing transistor circuit. 




FIG. 3. dc-characteristics of Ebers-Moll transistor model. 
Each characteristic curve corresponds to Vbb ranging from 1 
to 9 volts. Points along a particular curve corresponds to V cc 
ranging from to 60 volts. 

The work reported in this letter does not attempt 
to describe a method for constructing completely trans- 
portable behavioural models. A complete study of trans- 
portability would include an investigation of the models 
ability to describe a device operating at all bias points 
and subjected to all possible classes of drive signals. Be- 
cause, the drive signal a device "sees" is very much ap- 
plication dependent, and in lieu of knowledge of a spe- 
cific application, we will investigate a restricted notion 
of transportability. We demonstrate our methods for 
the device operating at many different dc-bias potentials, 
bu subjected to the same class of amplitude modulated 
(AM) drive signals. Thus, in this letter we investigate 



transportability across cd-bias potentials and avoid the 
issue of transportability across signal class. 
The class of AM drive signal we investigate is 

V = [(1 + Asm(ut))V c ] sin(nt), (1) 

with V c = 5V, uj = 5MHz, fl = 500MHz, and A = 4/5. 
This signal is added to the fixed dc-offset at Vbb ( see 
Fig. ||) with fixed dc-offset at V cc (see Fig. ||). The input- 
output data we work with is obtained by integrating the 
circuit equations for the Ebers-Moll model of Fig. [j] em- 
bedded in the amplifier shown in Fig. ^|. Approximately 
TV = 3000 input -output data points are generated from 
these equations by integrating them from to 3/zs, sam- 
pling every Ins. The inputs for the behavioural models 
are potential differences, (Vb, V ce ), the output data is the 
current, I c . In Fig. |we show examples of input-output 
data generated at the four different offset bias points in- 
dicated by the squares in Fig. [|. (We will refer to these 
data sets as pancakes.) 
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FIG. 4. Four input-output data sets (pancakes) along the 
characteristic curve corresponding to Vbb = 5V, indicated by 
the squares in Fig. ^[ 

II. NONLINEAR MODELLING 

The nonlinear behavioural models we reconstruct are 
radial basis function models. They attempt to predict 
the measured current, I c (the output) from the measured 
potential differences, Vb and V ce (the inputs). In passing, 
we remark that the inputs to the behavioural models are 
not the same as the AM signal used to drive the device. 

The models have the following form 

I c (t) = F[V b (t),V ce (t)}, (2) 
where F [u(t)] is a radial basis function with affine terms, 

K 

F[u{t)] =/3 + au(t)+5^w i #||c i -u(t)||). (3) 
i=i 



Here, u(t) — [Vb(t),V ce (t)}, while 0, a, and u>i are pa- 
rameters to be estimated. The Cj's are called centres and 
also need to be estimated. The function d> is a Gaussian 



(r) = exp 



r 
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where cr is another parameter which is kept fixed at a = 4 
(corresponding to half the standard deviation of the V ce 
data) throughout our investigations. 

For many cases it may be necessary to embed the volt- 
ages, Vb and V ce , using their past as well as present val- 
ues in the function approximation of Eq. (^|). This will 
certainly be the case for complicated drive signals. How- 
ever, for the class of AM driving signals we examined, 
the device produces signals which are simple enough that 
embedding can be avoided. 

When using radial basis functions an important issue 
is how many centres should be used, and where should 
they be located. An effective method of choosing appro- 
priate centres from a large set of candidates is the subset 
selection method of Judd and Mees [j^]. This method 
attempts to find the best selection of centres from a set 
of candidates. It selects the best set by evaluating a de- 
scription length criteria. We will use this method (sub- 
ject to the Schwarz Information Criterion) to determine 
(in some sense) an optimal model structure using data 
collected at one bias point. 



A. Constructing Simple Models 

To illustrate the effectiveness of using radial basis func- 
tions to model the input-output relationship between (Vb, 
V ce ) and I c , consider the data from the bias point at 
Vbb = 5V and V cc — 18V. In Fig. [| we show the input 
space for this data. (The full input-output space cor- 
responding to this data is the second pancake from the 
bottom of Fig. |.) 



FIG. 5. The input space of the data generated at 
(Vbb, Vcc) = (5, 18)V. The five circles are the locations of data 
points chosen as centres by the modelling method. 

The modeling method was given every fifth data point 
as a candidate centre. 



Ilipr^^iiiiiiiii 



700 900 1 1 00 

time units 

FIG. 6. A section of the response signal and the pre- 
dicted signal of the nonlinear model reconstructed at 
(Vbb, Vcc) = (5, 18)1/. We see that the fit is exceptional with 
the errors too difficult to see. 

The modeling method kept the linear terms in Eq. (|^) 
plus five radial basis terms. The data points selected 
as centres are indicated by the circles in Fig. [|. Fig- 
ure H shows both the true I c output and the output pre- 
dicted by the behavioural model. At this scale the two 
curves are indistinguishable, with RMS — 6.13 x 1CP 4 
or [RMS/STD(I C )} = 0.003. Clearly, for this data set a 
very good fit is possible with a small number of centres. 
We argue that because the modeling method chose non- 
linear terms the data is better described by a nonlinear 
model than a purely linear one. 
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B. Constructing Transportable Models 

Given the effectiveness of our method at constructing 
nonlinear models at one operating point, our goal now is 
to deliver a single nonlinear model which "works" over a 
wide range of operating points. 

There are two approaches which can be followed. The 
first is to develop "training" data which covers as wide 
a range of the devices input-output space as possible, 
and then construct a model using this data. The second 
approach is to generate "identical" data sets at many 
bias points, reconstruct models at each bias point, and 
form a global model by interpolating between the recon- 
structed models. We have attempted both approaches 
and have, so far, found the latter approach to be more 
fruitful. Therefore, it will be described below. 



In general, the relationship between the modeling pa- 
rameters and the bias points is 

P = 0(V bb ,V cc ) 
a = a(V bb ,V cc ) 
u>i = Ui (V bb ,Vcc) 

Ci = Ci (V bb , V cc ) ■ 

We will determine these functional relationships by con- 
structing models from input-output data corresponding 
to a few different bias points. Once this step is com- 
pleted, parameter values at bias points not explicitly in 
the "training" data can be estimated via interpolation 
or extrapolation. In this fashion, we can deliver a non- 
linear model that is transportable over a wide regime of 
bias points. The information needed for simulating the 
dynamics of a device is the operating bias point (from 
which we determine the parameters of the model) and 
the input signal sequence. 

Referring once more to Fig. ^ we notice that the pan- 
cakes appear to change smoothly with changes in the op- 
erating bias point. In Figj7] we show that the location of 
the centres shown in Fig. pTslowly vary as we move along 
the characteristic curve corresponding to V bb = 5V. We 
have found that this behaviour is true for all of the bias 
curves in Fig. ||. Also, the data shown in Fig. || was gen- 
erated in an "identical" fashion at each bias point. That 
is, the input-output data point with index 50 (say) in 
one pancake corresponds to the data point with index 50 
in any other pancake. Thus, if the number of centres is 
fixed, and if the centres in each pancake have the same 
indices, then the nonlinear models constructed with these 
centres will have the same structure at every bias point. 
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FIG. 7. The evolution of the five centres chosen by de- 
scription length along the characteristic curve. We observe 
that their locations smoothly change with changes in the bias 
points. 



When constructing a transportable model we keep the 
number of centres, and their indices, fixed. The exact 
location of a centre is given by curves like those in Fig. [?]. 
Therefore, only a least squares estimate of the weight 
parameters is necessary. Implicit in our approach is the 
assumption that a smooth change in the bias points will 
result in a smooth change in the response signal, which 
will manifest itself as a smooth change in the estimated 
model parameters. 

This is indeed the case as can be seen in Figs. ||. Fig- 
ures ||(a) and |](b) show two ways of obtaining an aver- 
age estimate of the evolution of the constant weight (3 
along the characteristic curve corrsponding to V bb = 5V 



j(d) are the results 
1(e) and gf) 



(V cc = 5 - 50V). Figures g(c) and 
for one of the linear weights ai, and Figs, 
show the same for a nonlinear weight, uj^. 

To understand how these estimates were obtained, 
recall that at each sample bias point we generated a 
N = 3000 point data set. At each bias point we used 
two methods to divide this large set into five different 
600 point data sets. In the first method, the first 600 
points of the large set forms the first small set, the next 
600 points forms the second small set, and so on. A (five 
centre) radial basis mode was constructed for each of the 
smaller data sets. The five estimates for the radial basis 
parameters are given by the dotted lines in Figs. ||(a), 
(c) and (e). The solid lines in the figures are the average 
values of the parameter estimates, while the starred lines 
are quadratic polynomial fits to the average values. 

In the second method we decimate the data by taking 
every 5 th point to obtain the 600 point data sets. The 
rest of the procedure for obtaining the model parameters 
shown in Figs. f§(b), (d) and (f) is the same as that just 
outlined. 

The point of the above is to emphasize that the param- 
eter estimates themselves have "error bars" and although 
the quadratic fits to the average are bad in places, they 
nevertheless describe the parameter evolution quite well 
as our final results show. 



C. New methodology 

Our investigation leads us to propose the following 
methodology for producing a transportable model for the 
electronic device. 

1. Sample the space of bias voltages, and generate 
input-output data sets, in an "identical" manner, 
at each bias point. 

2. Choose one of the data sets to reconstruct a model 
using the subset selection and description length 
methodology proposed by Judd and Mees Q . 

3. For each input-output data set, reconstruct a non- 
linear model (using centres whose indices were de- 
termined in step 2) and store the estimated param- 
eters together with their corresponding bias points. 



4. Fit a polynomial relationship between the esti- 
mated parameters and the bias points. 

Models constructed using this methodology should be 
transportable over a range of operating points. 

To use the behavioural model to predict the input- 
output dynamics of the device at one particular bias 
point, one estimates the values of the parameters in the 
model by interpolating the relationships between the pa- 
rameters and the bias point. Predicted responses to drive 
signals are generated by inserting the drive signal into the 
radial basis model. 

We have used the above methodology to produce a 
transportable radial basis model of the circuit discussed 
above. Behavioural models were obtained from data gen- 
erated at twelve bias points (Vbb = 3, 6 and 9, V cc = 5, 
20, 35, and 50) corresponding to four points on each of 
three characteristic curves. The locations of these points 
are shown by circles in Fig. ||. Because Fig. ^ is a two di- 
mensional projection of a three dimensional input-output 
space there are only nine circles indicating twelve bias 
points. Bias points (3,5), (6,5) and (9,5) are actually 
indicated by the same circle, and similary for (6, 20) and 
(9,20). 

At each bias point an AM driving signal was added as 
described above. The indices for the centres shown in 
Fig. H were used for all of the models, and a quadratic 
polynomial was used to describe the relationship between 
the model parameters and the bias points. 

The transportability of the models was tested using 
input-output data generated at nine bias points. These 
points are indicated by triangles in Fig. |3J. Once again 
and for the same reasons as above, bias points (3.5, 10), 
(5.5, 10) and (7.5, 10) are represented by the same trian- 
gle. (Similary for (5.5,25) and (7.5,25).) We emphasize 
that this test data represents bias points that are com- 
pletely different from those used to "train" the model. 
The results of these out of sample tests are shown in 
Table Q. We see that good performance is achieved and 
transportability of the model is clearly exhibited over this 
wide range of operating points. 



Test set 
(V bb ,V cc ) 


RMS 


RMS/std(7 c ) 


(3.5,10) 


0.0053 


0.0324 


(3.5,25) 


0.0047 


0.0143 


(3.5,40) 


0.0099 


0.0278 


(5.5,10) 


0.0052 


0.0412 


(5.5,25) 


0.0017 


0.0061 


(5.5,40) 


0.0129 


0.0361 


(7.5,10) 


0.0046 


0.0850 


(7.5,25) 


0.0050 


0.0296 


(7.5,40) 


0.0070 


0.0231 



TABLE I. The errors produced by the transportable model 
for twelve out of sample test sets. 



III. CONCLUSION 

We have introduced a new methodology for modelling 
the input-output dynamics of nonlinear electronic de- 
vices. Our study was particularly focused on producing 
a nonlinear model which was transportable over a wide 
range of operating points. We used radial basis func- 
tions to model the device. Transportability was achieved 
by finding a relationship between the models parameters 
and the operating points of the device. A key to the 
success was the fixing of model structure (i.e., fixing ap- 
propriate centres) so that interpolation of the parameters 
with respect to the bias points was possible. A topic for 
future work is to apply our methods to additional elec- 
tronic devices. A second is to investigate the "other half" 
of the transportability question, namely, can the response 
of the device subject to other classes of drive signals be 
predicted accurately by our models. 
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FIG. 8. The evolution of the estimated parameters along the characteristic curve corresponding to Vbb = 5V. The dotted 
lines show the estimates obtained from the chopped (a), (c) and (e), and the decimated (b), (d) and (f) data sets. The solid 
line is the average of the dotted lines. Also shown (starred) line) is the quadratic fit to the parameters using the bias points. 



